Realization of the purely spatial Einstein-Podolsky-Rosen paradox in full-field images 

of spontaneous parametric down conversion. 
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We demonstrate Einstein-Podolsky-Rosen (EPR) entanglement by detecting purely spatial quan- 
tum correla-tions in the near and far fields of spontaneous parametric down- conversion generated 
in a type-2 beta barium borate crystal. Full-field imaging is performed in the photon-counting 
regime with an electron-multiplying CCD camera. The data are used without any postselection, 
and we obtain a violation of Heisenberg inequalities with inferred quantities taking into account 
all the biphoton pairs in both the near and far fields by integration on the entire two-dimensional 
transverse planes. This ensures a rigorous demonstration of the EPR paradox in its original position- 
momentum form. 

PACS numbers: 03.65.Ud, 42.50.Dv, 42.50.Ar, 42.50.Lc 



I. INTRODUCTION 

In 1935, Einstein, Podolsky, and Rosen (EPR) pro- 
posed a Gedankenexperiment [1] involving two spatially 
separated but entangled particles. They showed that 
quantum mechanics predicts that these particles could 
have both perfectly correlated positions and momenta, 
in contradiction with the so-called local realism where 
two distant particles should be treated as two different 
sys- terns. Though the original intention of EPR was 
to show that quantum mechanics is not complete, the 
standard present view is that entangled particles do ex- 
perience nonlocal correlations [2, 3]. 

Recently, the realization and detection of entangled 
EPR states aroused much interest of the scientific com- 
munity fol-lowing a testable formulation of the EPR 
paradox introduced by Reid [4|, involving correlations 
of quadratures of twin beams. A recent review on the 
subject has been given in Ref. [5j. Furthermore, posi- 
tion and momentum entanglement of beams have been 
demonstrated at the EPR level by combin-ing squeezed 
light from two spatial modes, with measurements by ho- 
modyne detection in the temporal domain [H, [?J . 

Spontaneous parametric down- conversion (SPDC) pro- 
vides independent pairs of entangled photons that makes 
the systemvery close to that considered in the original 
EPR paper: the positions of photons 1 and 2 are detected 
in the near field, and their momenta correspond to the 
far field. Howell et al. havemeasured in both planes 
the probability distribution of the position of photon 2, 
conditioned by the detection of photon 1. The product 
of the conditional variances was found to be 25 times 
smaller than the limit for the product of variances for 
a single photon given by Heisenbergs uncertainty rela- 
tion. This impressive result was obtained by measuring 
tempo-ral coincidences between cross-polarized photons 
in type-2 SPDC. These photons were separated by a po- 
larizing beam splitter: for a fixed position of a narrow 
slit transmitting photon 1 to an avalanche photodiode, 



the level of coincidences was measured for each position 
of a similar slit transmitting photon 2 to a separate sim- 
ilar detector. Hence, in the words of Reid et al. [H, 
" detection events are only considered if two emitted pho- 
tons are simultaneously detected" . In this sense, they did 
not prospect the full EPR characteristic of SPDC but a 
monodimensional and point EPR paradox using postse- 
lected data. The same comment could be made on the 
recent realization by Leach et al. [9], where point detec- 
tion was used in one of the transverse dimensions of the 
field. 

We present here a full spatial demonstrationof 
EPRsteering by imaging highly spatially multimode 
type-2 SPDC with an electron-multiplying CCD (EM- 
CCD) camera. This ensures taking into account each 
detection event that occurs in a relatively long exposure 
time compared to the laser pulse duration and, more im- 
portantly, compared to the coincidence time detection 
used in the experiments that use correlation data detec- 
tion. Hence, even if there are losses and false-positive 
and -negative events, these spurious events are random, 
which is fundamentallydifferent fromconsideringaprior- 
ithat pairs are correlated and detecting only in the tem- 
poral and spatial gate where the twinphoton arrives. Sec- 
ond, Heisenberg inequalities concern a sole system, hence 
it is meaningless to test these inequalities if the near field 
and the far field do not correspond to this unique system. 
By treating full-field bidimensional images of photode- 
tection and measuring variances in two orthogonal direc- 
tions, we assure a perfect correspondence between the 
subsystems involved in the near field and in the far field 
[5j, in contrast to point or one-dimensional (ID) detec- 
tion. Moreover we demonstrate not only an EPR paradox 
in two spatial transverse dimensions but also one using an 
isotropic criterion involving the whole system. Note that 
the use of type-2 phase matching allows us to spatially 
separate the idler and signal photons to be close to the 
real conditions of an EPR test of local realism. Detec- 
tion is performed with an EMCCD camera. The ability 
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of EMC CD cameras to reach the photon-counting regime 
with a very high quantum efficiency [10|, LUL| makes them 
useful in quantum optics, and they have already been 
used by our group to characterize quantum correlations 
in the far field of type-1 [H, [H and type-2 Q SPDC. 
Note that previous work has been done by other groups 
to detect the quantum correlation in the far field [15|, [l6| , 
using intensified CCD (ICCD) cameras that have a lower 
quantum efficiency. 



II. THEORY 

SPDC induced by a wide monomode Gaussian pump is 
a strongly multimode beam: the extension of the down- 
converted beam in the near field (image plane) is iden- 
tical to that of the pump in the limit of low gain and 
for a sufficiently wide and thin crystal, while the far-field 
(Fourier plane) extension is limited by phase matching. 
The etendue of the beam, i.e., the product of its trans- 
verse surface by the solid angle it subtends or the number 
of transverse modes in appropriate units, corresponds to 
the two-photon Schmidt number [17j. 

The spatial extension of a mode in either the near or 
the far field is proportional to the inverse of the full beam 
extension in the other plane. For single-photon imag- 
ing, the laws of diffraction are equivalent to Heisenbergs 
uncertainty relation: a photon that can be localized in 
one mode of the near field, for example, by traversing 
an aperture of the size corresponding to the mode, will 
be detected at a random position in the entire far-field 
diffraction pattern. However, the laws of quantum me- 
chanics state that a pair of signal-idler photons will be 
detected either in the same mode in the near field or in 
opposite modes in the far field if no detection occurs in 
the other plane. Because the detection plane can be cho- 
sen at a time when causal interaction between photons is 
no longer possible, these correlations are not compatible 
with local realism, as demonstrated first in the EPR pa- 
per [l| , though compatible with Heisenbergs uncertainty 
relation since correlations cannot be measured in both 
planes for the same photon pair. 

We can describe the SPDC behavior as follows: for a 
detection of a photon 1 at i*i, the probability density of 
detection of a photon 2 at r 2 can be written as: 



N2 the corresponding quantity for photons 2, we have: 



p(r 2 \ r i) =p(r 2 ) + /(Ar) 



(1) 



where p(r 2 ) is the probability density of detecting a pho- 
ton of another pair (accidental coincidences) and /(Ar) 
is the probability density of detection of the twin photon, 
with Ar = ||r 2 ± ri||, plus sign indicating the far field 
(correlation of momenta on oppositemodes) and themi- 
nus sign the near field. Translational invariance, circular 
symmetry, and independence of the pairs (pure SPDC 
without further amplification) are assumed. Hence, if Ni 
is the number of photons 1 detected on a surface Si and 



N 2 )= f d 2 n f d 2 r 2 p(n,r 2 ) 



(2) 



= f d 2 n f d 2 r 2 {p(r 1 )p(r 2 )+p(r 1 )/(Ar)} 

J Si J S2 

Therefore, the probability of detection in S2 of the twin 
photon 2 of the photon 1 detected on Si is simply given 
by: 

F&)=f fr, /(Ar)= ^)-W^) (3 ) 
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If Si and S2 are the same size, this expression can be 
symmetrized and becomes the normalized intercorrela- 
tion function: 



F(S 2 ) = F(5i) 



(AW - (iVi) (N 2 ) 



(4) 



Themean in this equation can be estimated by spatial 
averages on the different pixels of the image for a fixed 
Ar, given by the intercorrelation of two "regions of in- 
terest" (ROIs) of an image, each one corresponding to 
one polarization of the SPDC. We will therefore obtain a 
nonlocal estimation involving all the light. Note that de- 
terministic spatial variations of the mean intensity do not 
preclude the validity of these spatial averages, inasmuch 
as the width of the intercorrelation function is smaller 
than the scale of this deterministic variation. Indeed, 
the covariance signal-idler for a region formed by inde- 
pendent area is the sum of the covariances of each area, 
just as the mean for the region is the sumof themeans for 
each area. Hence, if the ratio between the covariance and 
the mean intensity does not depend of this mean, it will 
be retrieved by spatial averaging even if the mean varies 
spatially. Because of the weak signal-to- noise ratio, we 
proceed to an additional statistical average on different 
images taken at different times for the same system con- 
figuration. 

For independent pairs, the quantity in equation (J2J) can 
be expressed as a function of the variance of the difference 
between Ni and N2: 

( Nl ) = (N 2 ) = (TV 2 ) - (N,) 2 

dm - iv 2 ) 2 ) 



F(S 2 ) = 1 



(Ni + N 2 ) 



(5) 



The physical quantities used to test the EPR viola- 
tion of Heisenberg inequalities are the spatial variances in 
each dimension A 2 x and A 2 y. To define a unique global 
criterion, we use the mean value of the two variances (i.e. 
half of the mean of the squared distance) : 



A 2 A 2 x,p x + A 2 y,p y 
A r, p — 



(6) 



Thus, by introducing the Heisenberg inequalities in the 
product of near-field and far-field variances, one gets 

A 2 rA 2 p = i(A 2 x + A 2 y)(A 2 Px + A 2 p,) > ^ (7) 
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We have used the fact that A 2 xA 2 p y + A 2 yA 2 p x > l\. 
Note that this two dimensional diffraction limit will only 
be reached by isotropic 2-D gaussians. 



III. EXPERIMENT 

The experimental setup is represented in Fig. [TJ The 
pump pulse at 355 nm provided by a passively Q-switched 
Nd:YAG laser (mean power: 27 mW, pulse duration: 
300 ps, repetition rate: 1 kHz), illuminates a 1-mm long 
type 2 BBO nonlinear crystal. The far-field image of the 
SPDC is formed on the EMCCD in the focal plane of a 
37 mm focal aspheric lens: Fig. QJa). In the near field 
configuration, Fig[]Jb), the signal and idler photons are 
separated by a Wollaston prism of 1.5 of angular sepa- 
ration positioned around the Fourier plane. The plane 
in the middle of the BBO crystal is imaged on the EM- 
CCD plane by a second identical aspheric lens, in order to 
minimize the walk-off effects [18|. The transversal mag- 
nification has been checked: 7 = 1.003 ± 0.005. The 
back-illuminated EMCCD camera from Andor Technol- 
ogy (model iXon+ DU897-ECS-BV) has a quantum effi- 
ciency greater than 90% in the visible range. The detec- 
tor area is formed by 512x512 pixels, with a pixel size of 
Spix = 16 x 16/im 2 (i.e. 0.46x0.46 mrad 2 in the far field 
after division by the focal length). We used a readout 
rate of 10 MHz at 14 bits and the camera was cooled to 
-85° C. Measurements were performed for a crystal ori- 
entation corresponding to non critical phase matching at 
degeneracy, i.e. collinear orientation of the signal and 
idler Poynting vectors in the crystal [Tjj. Photon pairs 
emitted around the degeneracy are selected by mean of 
a narrow-band interferential filter centered at 710 nm 
(AA=4 nm). As in [13| and [14[, the photon-counting 
regime is ensured by adjusting the exposure time in such 
a way that the mean fluence of SPDC was between 0.1 
and 0.2 photon per pixel, in order to minimize the false 
detections [Tcj. Moreover, the use of pump pulses with 
300 ps duration (much longer than the coherence time of 
SPDC) and an exposure time of the EMCCD of 10 ms 
(i.e. 10 laser shots) allow the excess noise to be limited 
by increasing the number of temporal modes [20] : the 
mean number of photons for one spatiotemporal mode is 
less than 10 -3 , in good agreement with the hypothesis of 
pure spontaneous parametric down conversion, without 
any stimulated amplification. 



A. Near field 

We show in Fig. [2] (a) the sum of the 10000 near 
field images. The two SPDC patterns are clearly visi- 
ble, with inhomogeneities and hot spots due to defaults 
on the crystal . The ROIs corresponding to either the sig- 
nal or the idler, large enough to encompass all the light 
for each polarization, have a common area. As a con- 
sequence, the inter correlation function exhibits a strong 




FIG. 1. Experimental setups, (a) Far field, (b) Near field. 
Brown arrows give the polarization directions 




FIG. 2. Intensity distributions, (a) Near field, (b) Far field. 
In both figures rectangles denote the ROIs used for calcula- 
tion. 

autocorrelation peak, as can be seen in Fig. [3j but the 
intercorrelation peak, due to quantum correlations, and 
the autocorrelation peak are clearly distinguishable. We 
use a Fourier algorithm without any zero padding to com- 
pute the intercorrelation, which is equivalent to a peri- 
odisation of the images. Figure |4] presents the near-field 
intercorrelation peak ob-tained from ROIs taken in the 
same image and a witness intercorrelation obtained from 
ROIs taken in two successive images. The absence of any 
non-negligible intercorrelation value in the second image 
shows that the inhomogeneities in the crystal do not cre- 
ate any deterministic intercorrelation pat-tern. The first 
image exhibits a weak intercorrelation vertical line (more 
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FIG. 3. Near- field complete intercorrelation in log scale. 
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FIG. 4. (a) Near- field intercorrelation function, (b) Corre- 
sponding witness intercorrelation. 



visible in [3] because of the log scale), which is due to the 
spreading of the autocorrelation peak induced by smear- 
ing in the gain register of the EMCCD camera. The 
existence of such an artifact again shows the importance 
of sufficiently separating the signal and idler photons on 
the camera to avoid any superposition of the autocorre- 
lation and intercorrelation peaks. By fitting the intercor- 
relation peakwith a two-dimensional Gaussian function, 
we estimate the inferred near-field standard deviations in 
pixels: 

Ax = 1.53±0.07, A?/ = 2.2±0.1, Ar = 1.89±0.09 (8) 

By integrating thefitted curve, we alsoobtain the total 
quantum correlation coefficient in the near field, R n = 
5 x 10 -2 . This coefficient corresponds to the intercor- 
relation of pixels much larger than the coherence area 



0. 



B. Far field 



Figure [2] (b) shows the sum of 10000 images in the 
far-field configuration. The intercorrelation function ob- 
tained in the far field is presented in Fig. [5] Note that the 
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FIG. 5. (a) Far-field intercorrelation function, (b) Corre- 
sponding witness intercorrelation. 



anisotropy of the peak is mainly due to the anisotropy of 
the shape of the pump. However, as predicted by simula- 
tions, an enlargement exists due to the imperfect degen- 
eracy of the photons wavelength. This enlargement is it- 
self anisotropic and, as predicted by the theory, is greater 
in the walk-off direction which separates the two fluores- 
cence spots in the far field (vertical direction on each 
image presented here). The experimental results are in 
agreement with this phenomenon since the enlargement 
of the intercorrelation peak is greater in the vertical di- 
mension than in the horizontal one [14[. We finally find 
the inferred standard deviations in far-field pixel units: 

Ap x = 2.35 ±0.08, Ap y = 1.85 ±0.07, Ap = 2.11±0.07 

(9) 

The total quantum correlation coefficient in the far field 
is R p = 4.4 x 10 -2 . In agreement with Eq. [5j we have 
experimentally shown in [14[ that (1 — R p ) is equal to the 
variance of the difference between areas greater than the 
coherence cell, expressed in shot noise units. 



IV. EPR VIOLATION OF THE HEISENBERG 
INEQUALITIES 

We are now able to test the violation of the Heisenberg 
inequalities by using the inferred quantitieswe havemea- 
sured. In the horizontal direction, we have: 



A 2 xA 2 p x 



1.53 x 2.35 27T ' Spix h 



/•A 

(0.048 ±0.008)ft 2 < 



giving a violation factor of 5.2 ± 0.8. 
In the vertical direction, we find: 

A 2 yA 2 p y = (0.06 ± 0.01)ft 2 < 
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(10) 



(11) 



giving a violation factor of 4 ± 1. 

And finally, by using the complete statistic of the flu- 
orescence, from ([6]) and ([7]) one obtain: 



A 2 rA 2 p = (0.06 ± 0.01)ft 2 < 



(12) 
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which gives a violation factor of 4 ± 1. 



V. CONCLUSION 

We have demonstrated a purely spatial EPR para- 
dox by using a full-field and direct detection method. 
The violation of the Heisenberg inequalities for the in- 
ferred quantities has been demonstrated for the whole 
system involving the two spatial transverse dimensions. 
By recording all the photon pairs generated by sponta- 
neous parametric down-conversion in the near field and 



the far field of the same system, we make fewer sup- 
plementary assumptions than previous papers that in- 
volve postselection or homodyne detection. Hence, this 
demonstration of one of the most fascinating phenomena 
of quantum mechanics is made in the form closest to its 
original formulation. 
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